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order electroweak and QCD phase transition has been discussed recently with the 
suggestion that detonations are stable. We examine here the case of a detonation 
more carefully. We find that in front of the bubble wall perturbations do not grow 
with time, but behind the wall modes exist which grow exponentially. We briefly 
discuss the possible meaning of this instability. 



PACS number(s): 98.80. Cq, 95.30.Lz 
*email: abney@rainbow.uchicago.edu 



A first order phase transition involves the nucleation of a bubble of one phase within 
a medium of the original phase. If the bubbles are large enough, they expand, collide and 
coalesce until the original phase has been completely replaced by the new phase. Though 
many of the details of this process are not definitively known — the expansion rate of the 
bubbles for instance — understanding the dynamics of such transitions may provide valuable 
insights into the conditions of the early universe. 

Of particular interest are the electroweak (EW) and QCD phase transitions and the 
possible ramifications towards the generation of a baryon asymmetry in the EW case and 
the concentration of baryons in the QCD case. If the EW phase transition is first order, 
the possibility of baryon asymmetry generation arises [0]. Furthermore, if the asymmetry is 
created by interaction with the bubble wall, it becomes important to understand in detail 
the shape and structure of the wall 0. A first order QCD phase transition, on the other 
hand, may result in baryon concentrating effects 0. Crucial to this is the concept of phase 
separation |Q. The generation of a baryon concentration may depend on whether, and how 
effectively, the two phases mix during the period of bubble expansion. The stability of 
the bubble wall and the possible existence of turbulence are thus important factors when 
considering the effects of these phase transitions. 

Recent investigations @, |], P] of the hydrodynamic stability of the bubble wall restrict 
their attention to the case of a deflagration front, namely a bubble wall which is propagat- 
ing at a speed slower than sound relative to the old phase. Recent estimates support the 
assumption that the bubble wall will propagate subsonically 0. Even if these estimates are 
correct, however, there may exist situations which allow for the supersonic propagation of 
the bubble, i.e., a detonation front. It was suggested by Kamionkowski and Freese that 
the existence of an instability in the deflagration front could result in the acceleration of the 
front until it becomes a detonation, a phenomenon which is observed in laboratories studying 
combustion ^J. The paper by Huet et al. ^ analyzed the dispersion relation obtained for 
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a perfect relativistic fluid and concluded that there could be no perturbations which grow 
in time in the detonation case. They, however, only examine the region in front of the ex- 
panding bubble and not the region behind it. In this paper we employ the standard linear 
stability analysis used by others |^, ^ [T^ to examine the possible existence of instabilities 
both in front of and behind the bubble wall in the detonation case. 

Detonations and deflagrations consist of two different phases separated by a transition 
region referred to as the bubble wall. A detonation front propogates at a velocity greater 
than the local speed of sound relative to the original phase. Behind the wall, the medium 
has a velocity equal to or less than the speed of sound in the second phase (see Steinhardt's 
article |ll|] for a discussion of relativistic detonation and shock waves). We choose to move 



into the frame of the moving wall and position it at x = 0. To the left of the wall {x < 0) is 
the original {e.g., quark) phase, while to the right (x > 0) is the new {e.g., hadron) phase. 
In region 1 (x < 0) the fluid has a positive velocity greater than that of sound {vi > c^i > 0). 
In region 2 (x > 0) the fluid has a positive velocity less than that of sound (0 < f 2 < 0^2) [PI - 
The situation is illustrated in Figure 1. Because of the different velocities, the behavior of 
perturbations and their effect on the bubble wall will differ. Since perturbations cannot 
travel faster than the speed of sound, we expect that in region one they should be "swept 
away" as the fluid passes through the wall. This is in agreement with references ^. The 
situation behind the wall is considerably different, there being no a priori reason why one 
should expect perturbations to decay with time. 

We first consider the behavior of a relativistic perfect fluid in the two separate regions. 
Using the metric g = (+, — , — , — ) the stress-energy of a perfect fluid is T'^^ = wu^u'^ — pg'^" 
where we have taken c = 1, with e the energy density, p the pressure, and w = e + p the 
enthalpy density. The equations of motion are: 

u^'df.p + clwd^u^ = 



2 



wu^d^Uy - dyp + UyU^d^p = 



where — (7, 7^) and 7=1/ y/1 — v^. We take perturbations in the velocity and pressure 
of the fluid about a constant solution: 

P=PQ + SP 

V = Vq + Sv 
Vq = Vox 

Sv — 5vxX + Svyy. 

Keeping terms up to first order in Sp and Sv, we can write the equations of motion as: 
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where (1) comes from the conservation of energy and (2)-(3) are the relativistic Euler's 
equations. It is more convenient to write this system as the matrix equation: 
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This system may be solved to obtain a solution of the form: 
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W{x,y,t) = F{x)e 
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where 
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The Qj are constants, and Rj are eigenvectors. The qj are found by solving the characteristic 
equation for (H). Doing this, we obtain the dispersion relation 



^{qvo + uj)'' -{q + vouj f - (1 - vl)k^ 







(5) 



The solutions for q are: 



l)voUJ ± Cs{l - vl] 



\ 



(6) 
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These, are the same as the equations obtained in reference equations (32), (34)- 

(36). It will be of future interest to note that in the case where Vq = we obtain two 
solutions, qi from above and 

02 = ^ -uj H . (8) 



Thus, the solution for perturbations in a perfect relativistic fluid is: 
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where 
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for j = 2, 3. 



\ w-f^{uj + Voqj) J 

We now wish to examine whether there exist any unstable modes of W [i.e., modes where 
Im uj > 0) which obey the boundary conditions W —>■ as x ±00. First, consider region 
one (x < 0). In order to satisfy the boundary condition as x ^ —00, we must require either 
Im qj > or aj = 0. However, with some algebra, we can show from (|^) that if Im u > 
then Im qj < and therefore aj = for j = 1, 2, 3. That is, perturbations which grow in time 
cannot exist in front of the bubble wall. This is the conclusion reached in reference fBll. Let 



us next consider what occurs in region two (x > 0). Here, the x +00 boundary condition 
requires either Im qj < or aj = 0. We find, if Im u; > then Im gi 2 < and Im > 0. 
Thus, behind the wall we need require only that 03 = 0. That is, since we can satisfy the 
X — > +00 boundary condition without requiring 01^2 = 0, the solution for the perturbations, 
equation (j^), is not identically zero, i.e., perturbations behind the bubble wall that grow in 
time may exist. 

The above work, however, is not enough to prove that instabilities do, in fact, occur. It 
is still necessary to impose the constraints of the boundary conditions across the bubble wall 
to determine whether the instabilities do not contradict the relevant conservation laws. Let 
us use notation where subscripts of one or two indicate quantities in region one (x < 0) or 
two {x > 0) respectively. Also, we assume that there exists a perturbation of the bubble 
shape of the form: 

A{y,t) = De-^(^*+'=^). 

We require 

• conservation of energy: 

Wi-flvi = ^272^2 

• conservation of momentum: 

^17^? +Pl= W2llvl +P2-T f — - ^ ) A 

• continuity of transverse velocity: 

dA dA 

Vly + Vl— = V2y + ^2-^- 

Recalling that = in region one, the above equations may then be linearized to obtain 
in region two: 
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where F-t = 1 ± 7|6'2w| and 02 = I + l/cs2- Then ([T0|)-(p!2[) give us the perturbations just 
interior to the bubble, i.e., as x — > 0^. Since we require as to be zero for x > we are 
left with oi, 02, and D as undetermined constants. By matching the general solution of 
perturbations for a; > 0, equation (^, with equations ([T^)-(|T^, we obtain three equations 
for these three unknowns. We will always have a solution for such a system provided that 
the determinant of the coefficient matrix vanishes. Let us write equations ([lO|)-(|l^) as 
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Matching solutions (|^) with ([I3D we get: 
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This equation may be rearranged into a single 3x3 matrix equation where the columns of 
the matrix are the vectors -Ri, R2 and Y . 
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By taking the determinant of the above 3x3 matrix we obtain an equation in u;, g2, a/w2, 
fi, f2, and Cs2- 
1 
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Combining equations (^ and ([T^) and solving for lo we can eliminate q2 and obtain an 
equation for as a function of fi, f2, 0^2, (ylw2 and k. By studying this equation we 



6 



can determine if there exist any cases where uj has a positive imaginary solution, thereby 
estabhshing that the detonation has a self-sustaining instability. 

Of particular relevance is the case of a Chapman- Jouget detonation. This occurs when 
the velocity behind the wall is equal to the speed of sound (f2 = 0^2) • Steinhardt showed that 
any detonation with spherical symmetry must be of the Chapman- Jouget type [O]. It has 



also been postulated that any "naturally" occuring detonation will meet this condition [|To| . 
In order to examine this case we use the value for given in equation (§). It is interesting 
to look at the case where the surface tension is zero (cr = 0); there are four solutions to cj, 
two of which are positive imaginary. The solutions are: 



uj = ±i\ - — ^ — Cs2k 

V 1 - Cs2 

The full equation with surface tension cannot be solved analytically. Instead we look at 
the behavior of u under different limits. In the long wavelength limit {k 0), keeping terms 
up to second order in k, equation (JT^) becomes: 



^4^(1 - ct)u;' + I'^^ict - 3cl + 2)u^ + ^^^(-ct + hcl- l))kW 

2W2Cs Vl W2Cs 2 

+ ii-clvi -cl + CsVi + 2)c^ ^^ = (15) 

Notice that the quantity a/w2 provides a natural length scale. In Figure 2, the behavior of 
Im uj was plotted as a function of k with values of f 1 = 1.5 Cs and Cs = 1/ V^. The particular 
choice for vi is arbitrary; the behaviour of Im uj is not altered significantly for different 
values. We have chosen a solution which matches with an unstable mode in the zero surface 
tension case, and we see that the rate of expansion increases slightly slower than linearly 
with the wave number, an effect of the surface tension. The discontinuity which appears as 
k increases is an artifact of the small k approximation; its location approaches the origin as 
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the surface tension is increased. That is, with a larger surface tension we must go to a larger 
wavelength in order to insure that the k —>■ approximation is valid, as we expect. 

Next, consider the short wavelength limit. As k goes to infinity only higher orders of k 
contribute and we get: 



c,A;^u; + i(fi -c,)c^A;^ = 0. (16) 



- C^)^ 4 5 4 



2W2 

Figure 3 shows Im c<j as a function of k with the same values as in Figure 2. Here we notice 
that as k gets larger Im uj quicky approaches a constant value of approximately 0.144w2/o"- 
Unlike in the deflagration case, where the surface tension results in there being a lower limt 
cutoff in the wavelength of instabilites, with a Chapman- Jouget detonation the growth rate 
does not drop to zero with shorter wavelengths, but rather reaches a maximum positive 
value. 

In order for these instabilities to be dynamically relevant the mode with the fastest growth 
time scale must be less than the time scale associated with the phase transition. The QCD 
transition lasts a time tu ~ 10"^ s and has a value of cr/w2 ~ 1 fm [Q. We showed above that 
Im UJ reached a maximum value of ~ 0.14^2/0", corresponding to a time of ~ 2.3 x 10~^^ s. 
Thus, there is ample time for the instabilities to mature. 

The scales associated with the EW phase transition are much smaller. In this case, 
using the formulas and parameter values as given in reference [r2|, we have a ~ O.Ogr^ and 
W2 ~ 40Tjf where is the critical temperature. With a critical temperature of 150 GeV 
we get cr/ti'2 ~ 3 X 10"^ fm. This leads to a maximum value of Im of ~ 5 x 10^ fm^^ 
corresponding to a time of ~ 6.7 x 10~^^ s. Since the phase transition lasts a time ~ O.OOStj^ 
|l^, where the Hubble time is tu ~ 10^^^ s, we see that there is sufficient time for the 



instabilities to grow. 

Our general picture of the instabilities, then, is the following. For very large wavelengths 



the instabilities grow only very slowly with time, vanishing as A approaches infinity. As the 
wavelength decreases Im uj rapidly increases with smaller A until a maximum rate is reached. 
Thus, unstable modes exist at all wavelengths. Even though these calculations were done 
for a Chapman- Jouget type of detonation, it is unlikely that the existence and time scale of 
the instabilities is so sensitive to the velocity behind the wall that this picture would be very 
much altered should V2 < Cg. Additional work beyond the scope of this report would need 
to be done, however, to verify this assertion. 

Note that these instabilities do not exist in front of the bubble wall, but rather behind it 
and at its surface. What, then, are the ramifications on the bubble wall and the fiuid inside 
it? To answer this, let us examine similar types of instabilites which exist in the laboratory. 
D'yakov and others have calculated the stability condition for shock waves in a classical 
non-relativistic fiuid. Though a shock wave and detonation are not identical, the fiuid 
dynamics are quite similar. It is possible to establish conditions in the laboratory which 
violate this stability requirement; that is, under certain circumstances there are unstable 
modes which exist on the surface of and behind the shock wave. Thompson et al. [TB| have 



recently done experiments which violate this shock stability condition. In their report they 
show a series of photographs where we see shock waves transform from planar to a billowy 
cloud-like surface as the stability condition is violated. The surface of the wave is travelling 
faster than the speed of sound, and the instabilities do not propogate forward from the shock. 
The transformation of the shock surface from planar to irregular is described by them as a 
"transition to turbulence" of the fiuid behind the shock. 

It is possible that a similar situation exists with the instabilites described above. The 
growth of perturbations behind the bubble wall may result in turbulence and a highly ir- 
regular surface. Such effects may entail an alteration in the picture of the phase transition 
proceeding through the growth of uniform spherical bubbles. This may very well be relevant 
when considering the possibility of baryon generation and concentration in the EW and QCD 
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phase transitions respectively. 

Whether detonations actually do arise, though, is not yet known. A more highly first 
order transition would likely result in an increased possibility of detonations. Another mech- 
anism could be the existence of instabilities in deflagrations as desribed in reference |^ . The 
above analysis, however, is limited to the linear regime and does not take into account any 
non-linear effects that may arise. The result of such effects may be to stablize the pertur- 
bations, as has been observed with deflagrations in the laboratory ||^. It remains to be seen 
whether non-linearities are indeed important. 

I would like to thank E. Kolb and G. Starkman for helpful discussions and comments. 
This work was supported in part by the University of Chicago Physics Department, the DOE 
and NASA grant NAGW-2381 at Fermilab. 
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Figure captions 

Figure 1: Schematic of the bubble wall in the frame of the wall. The wall is located 
at X — with the old phase {e.g., quarks) in the half-plane x < and the new phase {e.g., 
hadrons) in the half plane x > 0. 

Figure 2: Perturbation growth rate (Im a;) as a function of wave number {k) in the limit 
A; — > 0, plotted with parameter values Cg — l/-\/3, = 1.5 Cg in units of W2/cr. 

Figure 3: The same as in Figure 2, but in the hmit k ^ oo. 
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